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A simple way of calculating the transient response of a thin wire is presented. The method 
is based on (i) the asymptotic antenna theory for thin wires and (ii) an expansion of the induced 
current in terms of the so-called natural modes. An explicit series representation of the transient 
induced current in terms of the resonances of the wire is obtained and this representation is used 
to numerically evaluate the wire current using a desk calculator. A comparison is made between 
the results obtained using the asymptotic theory presented here and those obtained from a numerical 
solution of a space-time integral equation. Good agreement between the results is found except 
for early times. The method points to a fast way of estimating the transient induced currents 
on more complicated structures such as a wire model of an aircraft. 


1, INTRODUCTION 


So many papers have been written about scatter- 
ing and radiation from thin wires since Pocklington 
[1897] presented his before the Cambridge Philo- 
sophical Society that one more paper on this topic 
is hard to justify. Almost all work done in recent 
years involving transient scattering and radiation 
from thin wires has been based on numerical cal- 
culations ona computer [ Liu and Mei, 1973; Tesche, 
1973; Sayre and Harrington, 1972; Miller et al., 
1973]. These numerical calculations tend to become 
very complicated and there are many cases where 
a simplified analysis yields sufficiently accurate 
results. A simple and physically appealing way of 
calculating the radiation field from a pulsed antenna 
is presented by Franceschetti and Papas [1974]. 
The simple results so obtained are valid only in 
a time scale which is large compared to the transit 
time across the antenna, and thus fail to give the 
correct early-time response. However, the early- 
time response of a linear antenna has already been 
evaluated by Latham and Lee [1970] in their study 
of the transient properties of an infinitely long 
cylindrical antenna. 

In this paper, a technique based on the natural- 
mode method [Marin, 1973] is used to obtain 
approximate analytical expressions for the time 
history of the current induced on a thin wire when 
itis excited either by anincident step-function plane 
wave or by a slice generator whose output voltage 


Copyright © 1976 by the American Geophysical Union 


149 


is a step function in time. The advantage of the 
method lies in the fact that all expressions are of 
simple analytical form and can be evaluated using 
a desk calculator. The method is currently being 
used to obtain the currents induced on a thin wire 
when a charged particle moves near the wire. It 
is also used to study the external resonances of 
crossed wires (a simplified model of an aircraft) 
and two parallel wires (staggered and nonstaggered). 


2. MATHEMATICAL FORMULATION 


The induced surface current density j on a con- 
ducting body satisfies the integro-differential equa- 
tion 


Ff «j= —se,nxnx E'”, 


Gjds’ 


Ss 


2 ojmmanx(s?/e- VV» | (1) 


where 
G(r,r’) = (4n|r — r’|)~' exp(—s|r — r'|/c) 


is the free-space Green’s function, E'"* is the 
incident electric field, n is the outward unit normal 
to the surface S of the scattering body, and the 
time dependence exp (st) is understood. It is known 
that the singularities in the complex frequency (s) 
plane of the surface currents induced on a finite 
sized, perfectly conducting body can be identified 
either as the singularities of the incident field or 
as poles [ Marin, 1973] whose locations are uniquely 
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determined by the shape of the perfectly conducting 
scattering body. These properties of the induced 
current have been shown using the magnetic field 
integral equation and they imply that the inverse 
operator £~'(s) is a meromorphic function of s 
and that the locations of its poles (the natural 
frequencies) are given by all those s, for which 
the homogeneous equation 


L(s,)+j=0 (2) 


has a nontrivial solution j,,. 

It will now be shown that the operator £(s) 
is a symmetric operator when operating on functions 
that are tangential to S, i.e., 


L£1(s) = L(s) (3) 


where the ‘‘transposed’’ operator / T(s) is defined 
from the identity 


(£(s)+f,g) = (f,£7(s) -g) (4) 


In (4), f and g are tangential to S,i.e.,n-f=n-g 
= 0 on S and the scalar product is defined by 


(f,g) -| f(r) » g(r) dS 
Ss 


Using the definition (1) of Y together with (4) and 
some simple vector-algebraic manipulations one 
obtains: 


(Lf « fg) a) {cox Ee (s?/c? — VV +) 
Ss 
f Gcr,r’) f(r’) as’|| -g(r) dS 
Ss 


= -{ [G(r,r’) V -f(r') V - g(r) + s2/e? 
SxS 


- G(r,r') f(r’) + g(r)] dSdS’ (5) 


In view of (5) it is clear that (.% + f,g) = (f, 2 + g) 
showing that £ is a symmetric operator, i.e., 
= £T, 

The fact that “~-' is symmetric (an immediate 
consequence of the fact that Y is symmetric) and 
the assumption that £~' has only simple poles 
(which has been substantiated in all cases investi- 
gated numerically) make it possible to obtain the 
following forced solution of (1) by using the Mittag- 
Leffler theorem [see Marin, 1973]: 


B= ~-€98 DS) [65 — 5,)- (2 (8, )insdn) | (E™ 5) Jn 


+ &,(s) +E] —e,s8(s) BE (6) 


where £' = d.¥/ds, Z,,(s) is an operator-valued 
polynomial of sand % (s) is an entire operator-valued 
function of s. The quantity j,, is the current distribu- 
tion of a natural mode and it has the characteristics 
of a standing wave. 

Some comments are in order here concerning 
(6). First, the expression (6) is valid provided there 
is no degeneracy, i.e., for each s, there is only 
one linearly independent current distribution j, that 
satisfies (2). The case of degeneracy can easily be 
incorporated into the series representation (6) of 
the induced current in the same way as degeneracy 
was incorporated into the corresponding expression 
in Marin, [1973]. However, it is left out here since 
degeneracy does not seem to occur in the case 
of a thin wire. Second, the solution of (2) gives 
both the exterior and interior (cavity) resonances. 
On physical grounds it is also clear that s,, is purely 
imaginary for interior resonances and that Re{s,} 
< 0 for exterior resonances. The fact that Re{s, } 
# 0 for exterior resonances is an immediate conse- 
quence of the following two observations [see 
Weinstein, 1969a]: (i) the eigenvalues of the scat- 
tering operator have poles at the exterior reso- 
nances, and (ii) the scattering operator of a lossless 
object is a unitary operator for real frequencies 
(s purely imaginary). Also, the interior (exterior) 
resonances do not contribute to the sum (6) if the 
sources of the incident electromagnetic field are 
inside (outside) the scattering object. In fact, when 
the sources are outside the object it is shown in 
appendix A that the scalar product (E'"*,j,) van- 
ishes at the interior resonances showing that the 
residues of the interior modes vanish. This in turn 
implies that the sum (6) can be limited to include 
only the external resonances. 


3. RESONANCES OF A THIN WIRE 


The vector equation (1) can be reduced to a scalar 
equation for the total axial current I(z) when the 
scattering object is a wire of length | and radius 


a2 52 
( _ =) K(z — 2’) 1(z') dz’ = —se, EP(z), 


tt) 


0sz<!l i (7) 


where 


aa 
K(z- z') -| (8m? R)-' exp(—sR/c) do 


0 


and 
R2 = 4a? sin? (6/2) + (z— 2’)? 


Equation 7 can be solved asymptotically in terms 
of the ‘‘antenna parameter’? = 2 In(I/a) when 
| >> a. From this asymptotic solution the natural 
frequencies of a thin wire are found to be [ Oseen, 
1914]: 


s, = (c/) {ing — O-' [InQ| n|aP) -— CiQnz) 
+ iSiQnm)] + O(Q-2)},n = +1, +2,... (8) 


where ['(=1.781...) is the exponential of Euler’s 
constant, and Ci(x) and Si(x) are the cosine and 
sine integrals, respectively. Also, to the first ap- 
proximation, the current distributions of the natural 
modes are given by 


T(z) = 2maj,(z) = sin(naz/) +OQ-!), OSz<1 
(9) 


Expressions for the natural frequencies which are 
asymptotically correct up to and including order 
Q.~? can be found in Hallén [1930] where expres- 
sions for I,,(z) up to order .~' also are presented. 

To get some quantitative information about the 
accuracy of the asymptotic expansion (8) we have 
in Figure 1 graphed three different representations 
of the fundamental natural frequency of a straight 
thin wire, namely (i) the asymptotic form (8) correct 
up to order Q~' (labeled Ist order approx.), (ii) 
an asymptotic form correct up to order 0.~? (labeled 
2nd order approx.), and (ili) the numerical results 
obtained by Tesche [1973]. We note that for a/1 
= 0.01 (O = 9.2), the natural frequencies calculated 
from these different methods differ about 20% from 
each other. We also note that the second-order 
approximation gives too large a value of the damping 
constant |Re{s,}| whereas the first-order approxi- 
mation yields a somewhat too large value of Im{s,,}. 
Since the convergence of the asymptotic series is 
doubtful [Hallén, 1930], and judging from the 
results presented in Figure 1, it is questionable if 
the accuracy of the approximate solution would 
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numerical 
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Re{si/c} 


Fig. 1. The fundamental natural frequency s, of a thin wire. 

The natural frequencies for a/1 = 10-1, 10-5, 10-4, 10-3, 0.01, 

0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1 are indicated 
in the figure. 


be improved by including the 0~? order term when 
a/|= 0.01. For that reason and also for the reason 
that the calculation of the 0~? order term is rather 
involved for more complicated structures, we 
choose here to include only terms of order 07! 
in s, and terms of order 1(=°) in all other 
quantities. 

Some comments are in order here concerning 
the expansion parameter = 2 In(1/a) that is used 
in this paper. This parameter is not the only expan- 
sion parameter that has been employed in the 
asymptotic solution of cylindrical antenna prob- 
lems. The expansion parameters used by Wu and 
King [1965] and Schelkunoff and Friis [1966] 
depend on the frequency and they may lead to 
More accurate expansions of the induced current 
in certain cases. Since this approach is intended 
to be a simple (but not necessarily the most accurate) 
way of solving transient thin-wire problems, we 
find © to be the most convenient expansion param- 
eter. It should also be mentioned that expressions 
similar to (8) for the complex resonance frequencies 
of a thin wire have been calculated by Weinstein 
[1969b] who used the Wiener-Hopf technique to 
obtain approximate expressions for the induced 
current of a finite length wire. 
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4. CURRENT RESPONSE OF A THIN WIRE 


Using the approximate expressions (8) and (9) 
one can evaluate certain scalar products in (6) to 
get the following expression for the induced current 
on a thin wire: 


4ac Ss t 
—_—_—_— i "fl 
Zit 2 {— — S,) i ghar es) 


» EX (2',8) ac | sin (nwz/l) + O(0~*) 


I(z,s) = 


+ B&,(s)° ex} —e,s& (s) + Eire 


where the quantity E'"* (z,s) is the z component 
of the incident electric field. To determine the 
polynomial operator #,, and the entire operator ¥ 
we require that: (i) the expression (10) directly leads 
to the quasistatic solution for low frequencies, and 
(ii) in the limit of a very thin wire so that s, 
= inac/l, (10) is a solution of the differential 
equation 


(d?/dz? — s?/c?)I= -470-'segE™, O=72<1 


(11) 


with the boundary conditions I(0,s) = I(l,s) = 0. 
It is shown in appendix B that these two require- 
ments are satisfied when 2, = 0 and ¢ = 0. 

4.1. Gap excitation. When the wire is center- 
driven by a slice generator whose output voltage 
is a step function in time and has strength Vo, 
one obtains the following expression: 


(-1) 
2n+ 


1(z,t) = 


8V = 
= U(et - |z- 1/2) 5° 


sin [(2n 
Zo 1 [ 


n=0 


+ 1)7z/l] sin@,,,,¢) exp (—o,,,,, t) (12) 


where —o,, + iw, = S, and U(t) is the Heaviside 
unit step function. In Figure 2 we graph the time 
history of the current at z= [/4and z= 1/2 +6€,a 
<< e < |. Acomparison is also made with the results 
obtained by numerically solving a space-time do- 
main integral equation [Liu and Mei, 1973]. The 
approximate solution obtained here agrees with the 
numerical results obtained by Liu and Mei [1973] 


Asymptotic| Numerical 


ct/2 


Fig. 2. Step-function response of the current at z= 1/2 and z= 1/4 for a wire center-driven 
by a slice generator with output voltage V,. Also included for comparison are the corresponding 
results obtained by numerically solving a space-time domain integral equation. 


within 25% for ct/I> 1/4. It should be pointed 
out that the early-time response where the asymp- 
totic theory yields poor accuracy can be obtained 
from the results for an infinitely long cylindrical 
antenna excited at a 8-gap [Latham and Lee, 1970]. 

4.2 Plane wave excitation. When the wire is 
excited by a step-function plane wave whose direc- 
tion of propagation makes an angle 6 with the 
positive z axis and is so polarized that the electric 
field vector (strength E,) makes the angle 7/2 — 
6, 9 < 1/2, with the positive z axis, one gets the 
following asymptotic expression for the induced 
current: 


x 1 
I(z,t) = (8E,l/7OZ, sin 6) U(ct — zcos @) S | 
“in 


n=1 
NnwZ 
. ra {[sin@, t) — (-1)"sin@,t — nwcos6)] 


: exp(-0,0} (13) 


We note that the time origin is so chosen that the 
wavefront hits the wire end point z= 0 at t= 0. 


2Z,1 
E,2sin@ 
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The asymptotic expression (13) was used to numeri- 
cally calculate on a desk calculator the time history 
of the induced current at different positions on the 
wire -and at different angles of incidence of the 
plane wave. A comparison between these results 
and those obtained from a numerical solution of 
a space-time domain integral equation [ Liu and Mei, 
1973] is shown in Figure 3 for two angles of 
incidence, @ = 30° and 90°. It is observed in Figure 
3 that the asymptotic theory results exhibit faster 
oscillations than those of the numerical solution. 
The oscillations are due mainly to the fundamental 
resonance mode. An inspection of Figure 1 reveals 
that indeed, the fundamental natural frequency of 
the asymptotic theory has a larger imaginary part, 
implying faster oscillations than those obtained by 
numerically solving the integral equation. 


5. CONCLUSIONS 


Explicit expressions have been derived for the 
transient responses of a thin wire excited either 
by a plane wave or by a slice generator. The method 
has the great advantage of yielding explicit analytical 


3 4 


ct/2 


Fig. 3. Step-function response of the midpoint current for a wire illuminated by a plane wave 
with electric field strength E,. The cases @ = 30° and 90° are shown. Also included for comparison 
are the corresponding results obtained by numerically solving a space-time domain integral equation. 
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expressions for the induced current on the wire. 
This enables one to evaluate the current with ease 
at any point along the wire and at any time without 
having to resort to tedious numerical solutions of 
integral equations. The numerical results reported 
in this paper for the cases with the ‘‘antenna 
parameter’? Q = 10 are in good agreement with 
those obtained by numerically solving a space-time 
domain integral equation. Such an accuracy is quite 
adequate for many applications. The agreement has 
been found to be much better for thinner wires 
(i.e., © > 10). This is expected because of the 
asymptotic expansions used in the development of 
the theory. 

The time-domain expressions (12) and (13) are 
valid for an excitation waveform of a time step 
function only. It should be mentioned here that 
the method is also applicable to cases of other 
excitation waveforms. However, the excitation 
function E'"* (z,s) can have singularities other than 
a simple pole at s = 0 for more complicated wave- 
forms. Some care must then be exercised in evalu- 
ating the inverse Laplace transform of (10). For 
example, when the incident field is due to a charged 
particle that moves with a constant speed in a 
straight path parallel to the wire, then E'?s (z,s) 
has branch cuts along the positive and negative 
real semi-axes of the s plane. The results obtained 
for the ‘‘charged particle problem’ so far agree 
well with those calculated by the numerical solution 
of an integral equation, and they will be presented 
in a future publication. 


APPENDIX A 


In this appendix we show that the scalar product 
{E'"¢,j,,) vanishes at the interior resonance fre- 
quencies when all the sources of the incident field 
are located outside the scattering body. 

Let j,, denote the surface current density on S 
of an interior mode, and E,, H,, the corresponding 
electromagnetic field of this cavity mode in V, the 
region inside S. The incident field evaluated at s 
= s,, is denoted by Ei”*, H:"°. Some vector algebraic 
manipulations combined with Maxwell’s equations 
then give 


(E'™,j,) | Er. (nxH,) dS = -{ HT 
Ss RY 
i. E dV 
Vv 


+ (nxE,,) dS — i (Al) 


where i'"* denotes the sources of the incident field. 
The surface integral of the last expression in (A1) 
vanishes since the tangential electric field of the 
modes is zero on S. When all the sources of the 
incident field are outside S we therefore have 


(E35, 1,,= 0 


Equation A2 shows that the solution of the 
integro-differential equation (1) can be viewed as 
a superposition of an ‘‘interior solution’? and an 
‘‘exterior solution.’’ The interior (exterior) solution 
is due to sources inside (outside) the surface S. 
Another way of expressing the result (A2) is to 
say that all interior modes are orthogonal to any 
incident field whose sources are outside the scatter- 
ing object. It should also be pointed out here that 
the left-hand side of (1) is the same for both the 
exterior and the interior problem. It is therefore 
to be expected that the inverse operator £7 '(s) 
is singular at both the interior and exterior reso- 
nances of the body. The result (A2) shows, however, 
that there always exists a solution of (1) even at 
the interior resonances when the sources of the 
field are outside the object. The solution of (1) 
is of course not unique at these resonances although 
on physical grounds we define (6) to be the forced 
solution of (1) even at the interior resonances. 


(A2) 


APPENDIX B 


In this appendix we give two reasons for deleting 
the polynomial operators & ,, and the entire operator 
& in (10) for the case of a thin wire. 

To see that the expression (10) with 2, = € =0 
reduces to the quasistatic solution when s — 0 one 
lets s, = inac/l, n = +1, +2,... in (10) and obtains 


Sires f! 
I(z,8) = E™ (2', 8) 
ZA), 
ye cen ere Osineneyy) (BI) 
a s? + (nac/l)? 


The series in (B1), after dropping the s? term, can 
be summed as [Gradshteyn and Ryzhik, 1965]: 


» (nwc/1)~? sin (ntz’ /)) sin (nwz/1) 
n>0 


I—2)z'/2c%,z>2’ 
. z)z'/2c7,z>z és 


(l— 2')2/2c?, z< 2' 


Thus, (B1) becomes 


I(z,8) = 4tre, SQ! il z’ (1 — z/L) Ei (z',s) dz' 


t) 


I 
“| z(1 — z'/1) Bs sae | (B3) 


z 


and in the time domain we have 


0 Zz 
I(z,t) ~ 4meg 07! — il 2’ (1 — 2/1) Em (2',t) dz’ 
at LJ g 
Hy 

+{ z(1 — 2’/1) Eve wd | (B4) 
To show that (B4) is the correct quasistatic expres- 
sion for the induced current we note that the linear 
charge density +(z,t) induced on a thin wire by 
a quasi-electrostatic field with potential ® (x, y, z, t) 
is given by [Landau and Lifshitz, 1960]: 


7(z,t) = 47e,Q7' [@,(t) — ©(0,0,z,t)], O<z<l 
(B5) 


where ®,(f) is a function so that the wire has no 
net charge, i.e., 


1 
| t(z,Hdz=0 (B6) 


0 
Using the continuity equation we obtain the follow- 
ing expression for the induced current: 


Kn =-> | 7(z',t) dz’ (B7) 
at Jy 
which upon using the fact that Ei"(z,t) = — 
(0 /dz) ® (0,0,z, t) reduces to (B4) after integration 
by parts. 

To see that (10) with s, = innc/land 2, = @ 
= O satisfies the differential equation (11) we simply 
note that the Green’s function 


2 sin(naz/1)sin(naz’ /1) 


G(z,z')=— B8 
oe I > s? + (nac/ 1)? oe 

satisfies the differential equation 

(d?/dz? — s?/c?) G(z,2') = —8(z- 2’) (B9) 


and the boundary conditions G(0,z') = Gil, z’) 
= 0. In view of (B8) and (B9) it is now clear that 
the expression (10) with s, = innc/land 2, = & 
= O satisfies (11). It has to be emphasized, however, 
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that (11) together with the zero boundary conditions 
gives only undamped oscillations which are exactly 
the same as those given by (10) provided that 
radiation damping is neglected. 

Finally, we mention that the polynomial operators 
&,, and the entire operator # are not present in 
the series representation of the induced surface 
current density in the following cases: (i) all interior 
(cavity) scattering cases, (ii) scattering from a 
sphere, and (iii) scattering from a thin wire (in the 
asymptotic sense). Without a mathematical proof, 
we expect this to be the case for any finite sized, 
perfectly conducting scattering object. 
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